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ABSTRACT 

Obtaining reliable, matched photometry for galaxies imaged by different observatories represents 
a key challenge in the era of wide-field surveys spanning more than several hundred square degrees. 
Methods such as flux fitting, profile fitting, and PSF homogenization followed by matched-aperture 
photometry are all computationally expensive. We present an alternative solution called "synthetic 
aperture photometry" that exploits galaxy profile fits in one band to efficiently model the observed, 
PSF-convolved light profile in other bands and predict the fiux in arbitrarily sized apertures. Because 
aperture magnitudes are the most widely tabulated flux measurements in survey catalogs, produc- 
ing synthetic aperture magnitudes (SYNMAGs) enables very fast matched photometry at the catalog 
level, without reprocessing imaging data. We make our code public and apply it to obtain matched 
photometry between SDSS ugriz and UKIDSS YJHK imaging, recovering red-sequence colors and 
photometric redshifts with a scatter and accuracy as good as if not better than FWHM-homogcnized 
photometry from the GAMA Survey. Finally, we list some specific measurements that upcoming 
surveys could make available to facilitate and ease the use of synmags. 
Subject headings: methods: data analysis - techniques: photometric - galaxies: photometry 



I 1. INTRODUCTION 

Astronomy is entering an era of truly panoramic 
deep imaging surveys. Building on the legacy frorn 
■ the Sloan Digital Sky Survey (SDSS, lYork et alJ \2mi 
• lAbazaiian etall [2009I ). the coming years will see a 
rapidly expanding patchwork of overlappiiig imaging 
campaigns^^ with surveys like DE^, HSCFl KID^I, 
. VIKING^, and eventually Euclid and LSS'^B provid- 
ing hundreds of square degrees of new data on rapid 
timescales. Fast and reliable methods for combining 
information from these disparate surveys are clearly 
needed. 

A key challenge is measuring reliable matched photom- 
etry that samples, as closely as possible, the flux emitted 
from the same regions of a galaxy as observed across the 
many wavebands sampled by these data sets. For galaxy 
studies, the goal is to obtain accurate colors and spec- 
tral energy distributions (SEDs) over a large wavelength 
range in order to derive photometric redshifts (photo- 2;s) 
and physical properties like stellar mass (Af,) from tem- 
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plate fitting. But as extended sources, galaxies present 
particular challenges for matched photometry because 
their range in size, shape, and surface brightness is not 
easily accounted for by a single model or choice of aper- 
ture. 

Three basic methods have developed for obtaining 
matched photometry for extended sources from a set of 
multiband images with different seeing, scales, and other 
properties. The first and perhaps most common is PSF- 
homogenization followed by matched-aperture photome- 
try. Here, all of the overlapping images in a set of obser- 
vations are convolved to the PSF of the band with the 
worst seeing. The convolved images are then registered 
(often with interpolation to the same pixel grid) and the 
flux through a common aperture is measured on the pixel 
data in each band at the location of every source. This 
metho d has been studied in detail bv iHildebrandt et al.l 
()2012|) who stress the value of spatially-dependent con- 
volution kernels that are designed to make the PSF 
uniform across every image in all bands. The method 
has been successfully ap plied to a variety o f recent data 
sets including CO SMOS (|CaDak et al.l[200l and GAMA 
(|Hill et al.ll201lD . The disadvantages include the loss of 
information resulting from the PSF degradation, the in- 
troduction of correlated noise resulting from regridding 
and interpolation, and the fact that aperture flux mea- 
surements (especially for faint sources) are not ideal flux 
estimators. 

The second method for obtaining matched photome- 
try is "flux fitting." This approach is often adopted 
when the PSF sizes of the imaging data are strongly 
mismatched, for example when comparing photometry 
from HST to ground-based imaging or much lower res- 
olution near-IR data. In this case, the high-resolution 
image is first convolved to the lower resolution seeing. 
The pixel-by-pixel fiux of the "blurred" objects identi- 
fied in the high-resolution image are then compared di- 
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rectly to the pixel flux in the low-resolution image. A 
minimization scheme determines the flux ratio between 
the two bands for every high-resolution source, fitting all 
neighbors simultaneously to account for overlapping pro- 
files in the low-resolution imag e. This techn ique forms 
the basis of the TFIT pa ckage (iLaidler et al.1^ 0071 and 
the ConvPhot softw are (Ide Santis et al.ll2007i) and has 
been used widely (e.g-.fPapovich et al.ll2001l : lLabbe et al.l 
[20051: IGrazian et al.ll2006ir It has the advantage of mak- 
ing use of the available spatial information, but does 
assume — as all methods do to some extent — that the 
shape of the flux profiles are independent of wavelength. 
In this case, the shape is determined by the pixel-to-pixel 
flux in the single, high-resolution band. It can thus be 
subject to noise spikes and other anomalies. 

The final method is "profile fitting" which general- 
izes the "flux fitting" defined above. The most widely- 
used example is the photometry from SDSS, specifi- 
cally the ModelMa g photometry which is described in 
iStrauss et al.l (|2002[ ) , with more information available on 
the internelE3- Here a smooth intrinsic profile is deter- 
mined for every source, usually in the deepest band. In 
each of the other bands, this profile is convolved with 
the local PSF and scaled to the pixel flux in that band. 
Assuming the profile is wavelength-independent and cor- 
rect, this method maximizes the signal-to- noise (S/N) 
on the flux measurement and obtained colors. But the 
adopted proflle need not be a perfect match to the ob- 
served profile. In this case it can be thought of as a 
spatial weighting function applied t o the total fiux. This 
is the approach taken by iKuiikenI (|2008| ) for measuring 
the "Gaussian-aperture-and-PSF" (GaaP) flux, a seeing- 
independent flux estimator. Here the profile fitting is 
done to the shapelet basis function with fits obtained in 
every band. The shapelets can then be convolved analyt- 
ically to a common PSF under which the GaaP flux and 
resulting colors are defined. Since the shapelet fitting is 
performed in every band, however, the photometric scat- 
ter should be larger than in methods that fit all fluxes to 
a single intrinsic profile. 

Each of these methods has strengths and weaknesses, 
but all of them suffer from a major problem when con- 
sidering the size and pace at which new imaging data is 
becoming available. Multiband imaging data over thou- 
sands of square degrees can quickly sum to ~10 terabytes 
or more. If new data are continuously being added, keep- 
ing the matched photometry up-to-date can be extremely 
expensive computationally. In addition to the raw com- 
puting time, the human cost of mastering and manipu- 
lating so much raw data is significant. Instead, it would 
be valuable to have the option of measuring the photo- 
metric properties of sources in a given survey and filter 
one time only, but in a way that allows for a straightfor- 
ward comparison with other filters and surveys without 
revisiting the pixel data. 

In this work, we present an approach with this goal 
in mind that delivers "synthetic aperture magnitudes" 
(SYNMAGs), enabling very fast matched photometry with 
reasonable precision. Higher precision results can be ob- 
tained with a full analysis of all available pixel data (for 
example, using profile fitting methods), but SYNMAGs 

See http://www.sdss.org/dr7/aigorithms/photometry.html 
and http : //www . astro . Princeton . EDU/f hl/phot o-lit e . pdf 



provide a fast alternative until such an analysis can be 
undertaken. Our method works entirely at the catalog 
level. We require recorded parameters for profile fitting 
in at least one band, with measured aperture magnitudes 
(of arbitrary size) and PSFs recorded for the other bands. 
As a test case, we derive matched photometry from the 
combination of SDSS and the UKIRT Infrared Deep Sky 
Survey's Large Area Survey (UKIDSS/LAS) and show 
that its performance is as good, if not better, than PSF- 
homogenized methods. The SYNMAG technique is easily 
applied to any data set that overlaps with SDSS, but is 
applicable generally to surveys employing profile-fitting 
photometry in at least one band, such as CS82 (Erben et 
al., in preparation), HSC, and DES. We make the SYN- 
MAG code publicly available. In a foUowup paper (Higgs 
et al., in preparation), we present a second methodology 
specifically tuned to the SDSS-UKIDSS comparison that 
offers additional fiexibility. 

An overview of the SYNMAG method is presented in 
Section [21 with specific notes relevant to SDSS given in 
Section 12.21 The use of Multi-Gaussian Expansions is 
described in Section [231 while Section [3l presents various 
tests of SYNMAGs applied to real data. In Section [4] we 
conclude with some discussion of "application-oriented 
photometry" and describe how future surveys can pro- 
vide catalog information to facilitate the use of SYNMAGs. 
Appendix [X] provides recipes for applying the SYNMAG 
software to real data. 

2. METHOD 
2.1. Overview 

We begin with an overview of our approach. While 
our specific motivation is to apply this method to ob- 
tain matched photometry between SDSS and UKIDSS 
imaging data sets, the methodology will be appropriate 
for many other applications. As in nearly all treatments 
of multiband photometry including the use of matched- 
aperture photometry, we will assume there are no color 
gradients. This is obviously a poor assumption because it 
ignores how color depends on aperture size and obscures 
the meaning of an integrated color (see Section [4|), but 
is standard practice and simplifies the problem. 

Consider two imaging data sets that we will refer to as 
Da and Db- They account for two sets of filter bands, 
obtained in different conditions with different PSFs, per- 
haps on different telescopes. In the specific case we de- 
scribe later. Da will refer to the SDSS DR7 data set and 
Db to UKIDSS. Our method requires that profile fitting 
has been performed on the imaging data from at least 
one of the filter bands in one of the data sets. We will 
assume model profile fit parameters are available in the 
Da catalog and will choose a definitive band for these 
profiles, p, which is typically the deepest band in Da- In 
SDSS, p is defined as the r-band. Over all filter bands, 
i, available from both data sets, our goal is to obtain a 
reliable set of matched colors, Cpi, where Cpi = mp — rrii. 

We assume profile fits are not available in Db, other- 
wise matched colors could be obtained by comparing the 
profile fits in Da to Db- Instead, we assume that only 
aperture photometry has been measured and recorded in 
the object catalogs for Db- This is often the case be- 
cause aperture photometry is far easier to perform than 
profile fitting. 
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For a given band, i, we require a catalog measurement 
from Db of the spatially-dependent PSF, or PSFi. We 
then convolve the fitted 2D profile Ep from Da with PSF^ 
and integrate the result from R = Q to R = i?aper, where 
Rapev is the radius used to define the aperture photom- 
etry recorded in Z?s- Note that we have not assumed 
circular symmetry. While the apertures used are circu- 
laiFl. our method works with 2D profile fits that describe 
the varied, projected axis ratios of galaxies on the sky. 
We do assume that the true fiux is well described by 
smoothly declining profiles, however. In this way, we 
have "synthesized" the aperture flux that would have 
been measured in the p-band (from Da) under the same 
seeing conditions as the i-band (from Db) and within the 
same aperture as used to measure aperture photometry 
in the Db catalog. We use mp^sYN.i to refer to the re- 
sulting synthetic magnitude or "synmag." We can then 
construct the full set of colors across Da and Db as 

Cpi = Tip.SYN.i ~ ™i,apcr- 

In some cases, as in SDSS, colors internal to Da or 
Db may already be available and more reliable than the 
SYNMAG colors. These can be substituted for the syn- 
mag colors when available because the assumption of 
no color gradients means that, given two filter bands, a 
galaxy by definition has only one color. In practice, of 
course, synmag colors will differ from other techniques 
(e.g., ModelMag colors), but we emphasize that our goal 
is not to obtain the best set of matched colors, but the op- 
timal one given constraints on both computational time 
and the difficulty of downloading, understanding, and 
organizing large imaging data sets. Taking advantage of 
pre-computed colors internal to one of the data sets can 
lead to an overall improvement in the characterization of 
a galaxy's true colors. It is also possible to obtain colors 
from profile-fitting magnitudes derived in different bands 
and even with different profile shapes by comparing the 
flux integral of the fit profiles at a radius that encloses 
most of the light. If the assumed proflle shapes are sig- 
nificantly different, however, biases may be introducecF^ 
that depend on galaxy morphology and concentration. 

We can also use the set of colors, Cpi, to gener- 
ate a set of PSF-matched magnitudes. Let us assume 
that well-defined "total" magnitudes, mp^tot, are avail- 
able for the p-band in DA^ perhaps again via profile 
fitting. In SDSS, these could be defined as ^ModciMag 
or rcModciMag- Wc can then define a new set of total- 
magnitude, PSF-matched photometry for each band, i. 

Finally, we note that in the typical Db data set, aper- 
ture photometry is performed in several apertures of 
varying size. It is possible to compute SYNMAGs and 
resulting colors for each of these apertures. A weighted 
average could produce a better color measurement, and 
in the limit of very many apertures, this is equivalent to 
fitting the intrinsic profile, Sp, to the circularly-averaged 
radial flux profile in the band i. This moves closer to the 
construction of ModelMags in SDSS, in which Ep is fit- 
ted to the 2D pixel data in band i. We do not pursue 

In practice, our method could be extended to work with non- 
circular apertures. 

A sense of the magnitude of such biases is given by the dif- 
ference between Exponential and de Vaucouleurs magnitudes in 
SDSS, roughly 0.2 magnitudes ±0.1 for galaxies near r = 20. 



this extension here, however, in part because aperture 
photometry spanning only a few pixels near the centers 
of objects requires a careful treatment of the flux varia- 
tion within pixels. Even with well-resolved sources, these 
central aperture fluxes may not be computed reliably by 
many survey pipelines. Because they represent regions of 
high surface brightness, this may bias attempts to fit Ep 
directly to the tabulated aperture fiux measurements. 

2.2. Application to SDSS and UKIDSS 

Synthetic aperture photometry is valuable in a vari- 
ety of contexts, especially as an effective way to de- 
rive a first set of matched photometry ahead of a more 
detailed analysis involving the pixel data. Our origi- 
nal motivation, however, came from the desire to ob- 
tain optical through near-IR colors for galaxies targeted 
by the Baryon Oscfllation Spectroscopic Survey (BOSS). 
We desired matched colors between the SDSS imaging 
and UKIDSS. In this section we provide relevant details 
about these specific datasets, especially SDSS. These are 
not only important for understanding the ways we have 
tested our approach, they will be useful for applying syn- 
thetic aperture photometry to the overlap between SDSS 
and other data sets. 

The SDSS imaging is reduced and processed by 
the Photo pipeline, which has been packaged in 
sdssStools. The source detection, fiagging, profile fit- 
ting, and photometr y p erformed by Photo is described 
on the SDSS websit l"^ under "Algorithin s." The SDSS 
camera is described in'Gun n et al.l (|1998l ). the telescope 
system in iGunn et al. (200(|), and the filter system in 
lFukugita~et al. (1996). A technical summary is provided 
in lYork et al.. (2000) . Most relevant for our purposes is 
the profile fitting and definition of ModelMags, which 
are the default recommendation for obtaining extended- 
source colors within SDSS. 

A PSF-convolved, 2-dimensional de Vaucouleurs 
(deV), exponential (Exp), and PSF profile is fit to ev- 
ery detected source in every band, ugriz. For the deV 
and Exp profiles, the best-fit effective radius (along the 
major axis), axis ratio and position angle, and normal- 
ization in the form of a total magnitude — as well as the 
estimated errors on these quantities — are recorded in the 
PhotObj tables in the SDSS database. The standard 
form for the exponential profile, with r measured along 
the major axis, is given by 



^Exp(?') = /e.Expexp 



-'^E I - 1 

He 



(1) 



where ke — 1.67835, as needed to set the definition of 
Eg to the surface brightness at the half-light or effective 
radius. Re- In Photo v5_6 the Exp profile is truncated 
at large radii in the following manner: 



-^Exp 





^ixp'(0 - <! ^Exp[l - {r/Re - i?Y 3 < r/i?e < 4 



r/Re < 3 
3 < r/Re 
r/Re > 4 



(2) 

For the deV profile, "softening" is applied in the center 
so that the profile is given as 

For DR7, see ,http : //www ■ sdss ■ org/dr7/algorithms/index ■ html | 
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/dcv(r-) = /e,dovexp 



— 1 + 0.0004 



(3) 

with Kd — 7.66925. The truncation for the deV profile is 
accomphshed as follows: 



rSDSS 

-'doV 



(r) 




{r/Re - 7) 



212 



r/Re < 7 
7 < r/Re < 
r/Re > 



. (4) 

These definitions can be found in makeprof . c in the 
Photo package. We note that some examples of spe- 
cific SDSS profiles are contained in phFitobj .h, but the 
values given there for the total flux after accounting for 
the imposed truncation can only be reconciled if the trun- 
cation follows the form [1 — {r/Re — 3)^]^ " for Exp and 
[1 — {r/Re — 7)^]^ '^ for deV (that is, the terms are not 
squared). We instead take as correct the forms hard- 
coded in makeprof .c and described by Equations [2] and 
131 In terms of total magnitudes, the truncations lead to 
differences of 0.01-0.03 mag compared to the standard 
forms. Meanwhile, the applied profile softening at the 
center of the deV profile is more significant compared to 
the standard profile and can lead to differences of 0.05 
mag. 

The deepest band of imaging in SDSS is the r-band 
which is used to determine the profile-matched photom- 
etry known as ModelMags. These are defined as follows. 
For each object, the better of the r-band Exp or deV 
fits as measured in all the filters is taken as the default 
profile for that object. This profile is then scaled in nor- 
malization (but in no other parameters) to the pixel data 
in the other bands at the source location. The integral 
of the total flux of the resulting best fit in each band 
defines the ModelMag in that band. ModelMag photom- 
etry provides galaxy colors with less scatter than other 
estimators, a strong motivation for the SYNMAG approach 
which enables the use of ModelMags in conjunction with 
photometry matched to other imaging data sets. 

2.3. Mock Image Implementation 

In the next two sections we outline two methods for 
implementing synthetic aperture photometry. The first 
approach is more "brute force" in nature but provides 
a simple way of testing the more sophisticated Multi- 
Gaussian Expansion implementation that follows. 

For each source, we take the profile fit parameters de- 
scribing the intrinsic profile, Sp, from Da catalog. These 
parameters should include a characterization of the pro- 
file shape (e.g., exponential, de Vaucouleurs, or more 
generally a Sersic-n value), the effective radius, i?e, the 
profile normalization (e.g., the surface brightness, Sej at 
Re) and ideally the axis ratio. We then use these pa- 
rameters to create a mock postage stamp image of Ep 
with zero noise. The mock image is convolved with the 
spatially-dependent PSF of the target band in and 
the flux in the now convolved profile is summed over the 
relevant aperture(s) defined in Db^ returning the SYN- 
MAGs in the p-band for that source, mp^sYN,i- 

There are a number of obvious limitations to this im- 



plementation. First, our goal has been to achieve a very 
fast way of producing matched photometry. The image 
convolution here is the slowest step, and while we have 
avoided returning to the raw pixel data, the mock image 
implementation is not that different than working with 
the real data since both require convolutions performed 
on 2D images. It is not clear that the ease of this ap- 
proach justifies giving up the advantages of a full-fledged 
re-analysis of the imaging in Da and Db- 

A second problem compounds the issue. Because our 
starting point is the intrinsic profile, Sp, the pixel scale 
"resolution" of the mock image has to be high enough 
to accommodate the steeply rising interiors of profiles 
with small i?e, especially deV profiles or those with high 
Sersic-n values. In our brute force approach, we found 
the needed resolution simply by decreasing the pixel scale 
of the mock images until there were no differences in 
the SYNMAGs measured at different resolutions. In tests 
of massive galaxy profiles at z ~ 0.5 (see Section 13. ip 
we found convergence at a mock image resolution of 
roughly 0.01 arcsec per pixel. This resolution results 
in prohibitively large PSF convolution times that were 
the whole motivation of synthetic photometry. A more 
clever approach with an adaptive resolution grid would 
address this problem, but given the other limitations, we 
adopt a completely different and more elegant approach 
in the following section. 

2.4. Multi- Gaussian Expansion Implementation 

Synthetic aperture photometry becomes attractive 
compared to a reanalysis of imaging pixel data once it 
is significantly faster. Because the slowest step in meth- 
ods that work with pixel data is often PSF convolu- 
tion, we develop a method here which decomposes intrin- 
sic profiles and PSFs into multiple Gaussian expansions 
(MGEs). The fact that convolutions between Gaussians 
can be calculated analytically makes synthetic photom- 
etry with MGEs extremely fast, as long as the MGE 
decomposition of the profile and PSF need only be per- 
formed once for the full data set. We will show this is 
commonly the case. 

The power of MGEs in the context of smoothly de- 
creasing profi les common in astronomy was recognized 
early on fe.g., lTrumpler fc WeaveH[i953[ ) and applied rig- 
orously to flux profile s such as the King model and Mof- 
fat PSF by EendEsIIl ([i99l|). Th eir utility and a modern 
implementation was presented in lErnselleni et al.l ()1994[ ) 
and also discussed in lConnollv et al.l (|2000t ). 

2.4.1. MGE Decomposition 

Our flrst step is to describe the intrinsic profiles, Ep, 
that are fit to the profile-band data in the dataset. Da, 
by an MGE, that is, a sum of j Gaussian profiles: 



(5) 



where Aj is the normalization (central fiux density) and 
the width of each Gaussian is given by <Tj . The summed 
flux of each Gaussian is Lj — 2TTAja'j. As we describe 
below, it is possible to determine an MGE for a pro- 
file of arbitrary size and shape, a process that takes less 
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Fig. 1. — A sample Exponential (left) and de Vaucouleurs (right) profile, illustrating the MGE decomposition. The standard forms are 
shown as solid purple lines, while the SDSS-modified forms are shown in dark grey. Surface brightness units are arbitrary and scaled so 
that E(iJe) = Ee = 1.0. The MGE components for the modified profiles are overplotted as thin, light grey lines. They sum to give the 
dashed blue curves. 

than a minute. However, as long as the chosen intrinsic 
form of the profile scales with an effective radius (Re) 
and normalization (Ee), the (Jj and Aj of a single MGE 
determined for that profile form can also be scaled in the 
same way. For a defined profile shape, this means a single 
MGE can be scaled to match profiles with arbitrary size 
and flux. Determining the MGE for the intrinsic profile 
of every source in Da is not required, although a separate 
MGE is needed for each profile shape that was fit for. In 
SDSS, this means an MGE must be determined for both 
the Exp and deV profiles. For Sersic profiles generally, an 
MGE is required for each Sersic-rt value. From the point 
of view of SYNMAGs, it is helpful if the number of profile 
shapes fitted for in Da (i.e., the number Sersic-n values 
fitted for) is restricted to a modest number. The MGE 
descriptions of these shapes can be determined once and 
then used for any subsequent measurement of SYNMAGs 
involving the datasets in question. 

So far we have considered axisymmetric profiles defined 
in one dimension. Two dimensional profiles are easily 
described by a sum of 2D Gaussians. We let the profile's 
major axis be the x-axis and define the axis ratio, q = 
b/a, where a and b are characteristic sizes along the major 
and minor axes respectively. The MGE becomes, 
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Note. 

S(Re) = 



The components are normalized so that the sum yields 
= 1.0 in arbitrary units of surface flux density. The units 



Sp(a;,y) = ^Aj-exp 



2ai 



(6) 



after substituting ay = qox- For a given profile shape, 
any 2D representation of arbitrary size, flux, or axis ratio 
can thus be described by a single ID MGE. 

The MGE must now be carefully determined with an 
eye towards minimizing the number of Gaussian terms 
while preserving precision in the MGE approximation of 
the intrinsic profile. We find that a least-squares fit to 
the full (circular) 2D profile provides the best approxima- 
tion to the true surface brightness profile because the fit 
is then weighted by the intensity. MGE approximations 



for standard Exp and deV profiles as well as the modified 
profiles used in SDSS are given in Table [1] and their ID 
radial profiles are shown in Figure [T] In terms of enclosed 
flux, the number of terms used in these MGEs provide 
accuracies at the 0.01 mag level out to 10i?e and so we 
deem them adequate approximations for our purposes. 
More details on the optimal determination of these flts 
as well as more accurate representations employing ad- 
dition al Gaussian terms are provided in iHogg fc Lang! 

(HoTp). 

It is generally difficult to fit the imposed truncation in 
the modified SDSS profiles with Gaussians, although the 
central smoothing in the modified de Vaucouleurs profiles 
makes an MGE solution easier in the central regions. The 
MGEs begin to deviate from the true SDSS profiles at 
radii larger than roughly 3.5i?e for the Exp profile and 
6i?e for the Dev profile (see Figure [T|) . This introduces 
aperture fiux discrepancies of roughly 0.01-0.02 mag, and 
so we consider these MGEs adequate for our purposes. 

2.4.2. PSF Convolution with MGEs 
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The observed galaxy profile, Ip is the convolution of Ep 
with the PSF profile, i?psF- The Fourier transform, J^, 
of this convolution can be expressed by the product of 
Fourier transforms: 



T{^p ® i?PSF) = ^(Sp) J-(i?PSF) 



(7) 



The power of the MGE approach lies in the fact that 
the Fourier transform of a Gaussian is also a Gaussian 
and easily determined analytically. Fourier transforms 
are also linear, so that for each Gaussian term, j, 



^(I]p(r)) = ^ J A,exp|^- 



2T:ikr dr. (8) 



Referring to Equation[6l T(Yip{x,y)) can be expressed 



as 



J"(Ep(a;,y)) 



27r Aj cr| q exp 

j 



(9) 

For the derivations here, we assume that -Rpsp can also 
be decomposed into a series of circular Gaussian terms, 
fc, with coefficients Rk and standard deviations, <TB,,k- A 
2-component Gaussian PSF, often used to account for 
power in the wings of the PSF, is an obvious example. 
Taking the inverse Fourier transform of Equation [71 we 
can express the PSF-convolved "observed" profile, Ip, as 




where the coefficients Cj^k are given by 




' R,k 



'R.k 



(11) 



While we have assumed circular PSF components to 
this point, the SYNMAG approach and associated software 
handles non-circular, spatially offset, and even negative 
components, providing a flexible way to describe more 
comphcated PSF shapes via i?psF. 

2.4.3. MGE Aperture Flux and profile 

Equation 1101 provides an analytic approximation, Ip, 
of a galaxy light profile with an intrinsic shape measured 
in the p-band, but observed with the PSF of another 
band. In our case, the second band comes from a differ- 
ent dataset, Db for which only aperture photometry is 
available. To compare the photometry recorded in 1?^ 
to the aperture flux, Fp, that would have been measured 
in the p-band with the same PSF, we must integrate 
Equation [TOl over the aperture(s) used in Db- Because 
Ip is generally elliptical while the apertures are circular, 
a numerical integration is required. Working in Carte- 
sian coordinates, the 2D integral of Equation [10] over a 
circular aperture of radius R can be expressed as 



FAR) = 

EE / / 4C7,,fcexp 



2ct? 



exp 



-^jdxdy, 
(12) 



where c(x) — \/R^ 



and a represents the effective 



standard deviation of a PSF-convolved MGE term, such 
that al = {alj + a\ f?i, and al = {q^crlj + a^ k)- The 
second integral can be expressed in terms of the error 
function, so that a fast ID implementation of the nu- 
merical integral can be written as 



Fp{R) 



'iCj^kCTy^j - 



exp 




da;. 



(13) 



We can also use the MGE description of Ip to estimate 
the effective radius and axis ratio of the observed fiux 
profile. Working along the major axis — or equivalently 
assuming the axis ratio is 1.0 — we can define an effective 
half-light radius, i?e,inaj, analogous to the Re defined for 
the intrinsic profile, Sp. The integral of the 2D flux is 
easily calculated: 



L{x) 



j.k 



2TTCj^k(^x 



1 — exp 



X 



(14) 



Setting this equal to the (major-axis) "half-flux" given 
by ii/2,maj = J2''^^j:kO'x ^-ud solviug for X, one obtains 



R 



e,maj • 



In practice, this definition of half-light radius is 
not as intuitive as in the case of Sp because the axis ratio 
depends on radius (as we show below). Thus, the minor- 
axis half-light radius is not simply related to i?e,maj • We 
can alternatively define a second measure of the half-light 
radius, i?e,circ, that corresponds to the size of a circular 
aperture which contains half of the total fiux. Here, we 
integrate Eg nation [T51 until the enclosed 2D flux is equal 
to the true half-flux as given by L1/2 = J^'^^^'j-k^xO'y 

Turning to the effective axis ratio, while the intrinsic 
axis ratio is flxed, the observed value will be a function 
of radius. Seeing will round out profiles at a scale cor- 
responding to the PSF FWHM, but will have a smaller 
effect on elongated (or inclined) profiles at larger radii. 
The convolved axis ratio of each Gaussian term is given 

by 



qj,k - 



'R,k 



(15) 



'R,k 



To derive an estimate of the observed profile's axis ra- 
tio, (g), we can weight each term by approximately the 
half-flux: 



(a) 



(16) 
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Fig. 2. — Comparison between r-band FiberMags and SYNMAGs. Sources best fit by exponential profiles are compared in the left-hand 
column, while de Vaucouleurs fits are shown in the right-hand column. From top to bottom, the panels investigate magnitude differences 
as a function of intrinsic axis ratio (g), effective radius (Re), and FracDev, an estimator for how well the profile shape can be fit by a de 
Vaucouleurs profile. Residuals at the 0.01-0.02 mag level are probably caused by an imperfect treatment of the PSF. 



where the weights are given by Wj^k = Cj^kd'xO'y. 

3. TESTS WITH GALAXY PHOTOMETRY 

The motivation for SYNMAGs is fast matched photom- 
etry for overlapping data sets that span large areas of 
the sky and thus make a reanalysis of pixel data difficult 
and time consuming. Since we cannot use the stellar 
locus — a common way of testing photometry methods — 
we perform tests of the SYNMAG approach using galaxy 
magnitudes and colors. We first show that synthetic 
magnitudes are consistent with SDSS FiberMags. We 
then compare galaxy colors from synthetic photometry 
to "FWHM-homogenized" matched-aperture colors from 
the GAMA Survey. We consider the derived scatter in 
the color of the red sequence and in the photometric 
redshifts determined from these two sets of photome- 
try, finding that SYNMAGs perform as well or better than 
the standard but more laborious FWHM-homogenization 
method. The latter two tests require redshift information 
which we take from the Baryon Oscillation Spectroscopic 
Survey (BOSS). 

3.1. SDSS FiberMags 

In all five bands, and for every detected source, the 
SDSS Photo pipeline provides a "FiberMag" estimate 
of the 3" diameter aperture magnitude, assuming seeing 
of 2" FWHM. This is meant to approximate the flux 
measured by a spectroscopic fiber in SDSS-I/II. Since 
SYNMAGs can provide synthetic aperture photometry for 
an aperture of arbitrary size and with arbitrary seeing, a 
natural first test is to compare SYNMAG approximations 
of FiberMags to their actual values. 

With an eye towards future applications of SYNMAG 
photometry for BOSS galaxies, we perform this exercise 
usin g a sub- region of the SDSS Stripe 82 Coadd cata- 
log (|Annis et alJ [20T1I) with 355° < RA < 360°. We 



select extended sources with non-stellar "psfmag" col- 
ors and detections in the UKIDSS LAS catalog (DR4). 
This effectively selects galaxies with A/* > lO^^Af0 and 
0.2 < z < 0.7. We apply cuts on both SDSS and UKIDSS 
flags to ensure a clean sample and remove blended ob- 
jects by excluding SDSS sources with the blend flag set. 
The exact nature of the flag cuts do not influence our 
results. 

We have found that a precise understanding of the PSF 
is important for resolving discrepancies at the 0.05 mag 
level. In the SDSS imaging, the PSF can be approxi- 
mated by a 2-component Gaussian. The second compo- 
nent has a width that is twice that of the first compo- 
nent, but accounts for only 10% of the PSF power. In the 
Photo pipeline (see measureObj . c), the FiberMag is de- 
termined by convolving a sub-image by a single Gaussian 
kernel with cr^ = cr| g — cr^o^, where a2 is the value of 
a for a Gaussian with FWHM = 2.0" (i.e., 2.0/2.3548) 
and (Tmom is the adaptive second moment size of the lo- 
cal PSF (referred to as M_rr_cc). Thus, the FiberMags 
are not measured on images with a FWHM 2.0" Gaus- 
sian PSF but with a more complicated PSF profile with 
significant power in the wings. 

For the consistency check here, we do not attempt to 
recover the applied kernel or the original (and likely spa- 
tially varying) PSF of the underlying image. Instead, we 
approximate the final, effective PSF for FiberMags as 
a 2-component Gaussian (as described above), with the 
FWHM of the first term set to 2.05". Experimentation 
revealed that our uncertainty in this choice introduces 
systematic offsets at the 0.01-0.02 mag level. We then 
divide the sample into sources that are best fit by ex- 
ponential profiles versus de Vaucouleurs profiles (this is 
easily determined by whether the ModelMag matches the 
ExpMag or deVMag) to see how the SYNMAG MGEs and 
PSF convolutions perform in either case. 
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Fig. 3. — Restframe (r — J) color-magnitude diagrams of (mostly) red-sequence CMASS galaxies from the BOSS survey. The observed 
scatter for the GAMA FWHM-homogenized matched-aperture photometry is larger than for the SYNMAG approach described here. Contours 
indicate the location of much less luminous COSMOS galaxies in the same redshift range. 

The results are shown in Figure [2] which compares 
r-band magnitude differences between SYNMAGs and 
FiberMags as a function of axis ratio, size, and FracDev, 
a quantity that expresses how well a source can be de- 
scribed by a de Vaucouleurs profile. The SYNMAGs per- 
form well with small systematic residuals at the 0.01- 
0.02 mag level that probably result from our imperfect 
treatment of the PSF. We expect additional systematics 
from the discrepancy between the MGE profile approx- 
imations and the truncated SDSS profiles and from the 
fact that a single-component profile is not a good descrip- 
tion of most galaxies. For our purposes, these system- 
atics are less than the statistical uncertainties in most 
applications of galaxy photometry as well as the photo- 
metric calibration uncertainties inherent in comparisons 
of SDSS to other data sets. Finally, we note a small de- 
gree of scatter towards brighter FiberMags. A large frac- 
tion (^80%) of these objects were removed by restricting 
to non-blended sources. We therefore surmise that the 
remainder are objects with nearby neighbors that con- 
tribute flux to the FiberMag aperture or galaxies with 
bright sub-components that are not accounted for in the 
profile fitting. 



As des cribed in Driver et al.l ()2011l ) and further de- 
tailed by iHill et al.l (|2011h . the GAMA DRl photomet- 
ric catalog was constructed by first downloading the re- 
duced images in all 9 available bands — ugriz from SDSS, 
and YJHK from LAS — and convolving all imaging to 
a common seeing FWHM of 2.0". We refer to this 
as "FWHM-homogenization," and distinguish it from 
"PSF-homogcnization" in which the FWHM and de- 
tailed sha pe of the PSF are made identical in every band 
(see Hilde brandt et an [2012:V All images are registered 
to a common astrometry, and source detection is per- 
formed in the r-band, which defines the Kron aperture 
that is applied to the flux measurements in the other 
bands. In the G09 and G15 fields, the catalog is limited 
to a Petrosian magnitude of Tpct < 19.4, while in G12 



the limit is 



pet 



< 19.8. Further details and various tests 



3.2. Red Sequence 

Having verified the ability of SYNMAGs to reproduce 
aperture photometry measured in the SDSS r-band, we 
now investigate how well SYNMAGs provide matched pho- 
tometry across differing data sets. Our focus will be 
on regions of overlap between SDSS imaging and the 
UKIDSS/LAS. These data sets have been matched pre- 
viously using standard FWHM-homogenization with re- 
measured aperture photorn etry in two recent papers. 
iMatsuoka fc Kawaral ()2010f ) carry out this exercise for 
40 de g^ in Stripe 82. The GAMA survey priver et al.l 
120111) does so in three equatorial fields, totaling 144 deg^ 
(|Hill et al.ll2011[ ). Because the GAMA DRl catalogs are 
larger and publicly availableP^. we will use them to test 
the synthetic photometry method developed here. 

|http : //www ■ gama-survey ■ org/dr 1/YRlpubllc ■ php | 



of the performance of GAMA photometry are provided 
in IHill et all (I20TI . 

Our first test measures the scatter in the restframe 
(r — J) colors of so-called CMASS galaxies from the 
Baryon O scillation Spectroscopic Survey (BOSS, see 
lEisensteiiT e t al. 2011i), most of which lie on the red se- 
qu ence. The BOS S targe t selection is briefiy described 
in I Anderson et al.l ()2012D and will be further detailed 
in Padmanabhan et al. (in preparation). The CMASS 
selection is designed to target very luminous galaxies 
with z > 0.4 and > lO^-^M©. Most have red col- 
ors, although 20-30% show evidence for disk or irreg- 
ular components, and 15% have s omewhat bluer color s 
than a standard red-blue color cut ([Masters et al.|[20T]] ). 
We select 1203 CMASS galaxies in the three GAMA re- 
gions with secure spectroscopic redshifts. We apply K- 
corrections to infer restf rame absolute magnitudes using 
K-correct (jBlanton fc R oweis 2007) and apply no evo- 
lutionary corrections. 

SYNMAGs are derived using the met hod described in 
previ ous sections from the SDSS DR7 (jAbazaiian et al.l 
12009( 1 and UKIDSS/LAS publicly-available photometric 
catalog^. Detections are required in all bands. We 

Special care is required in querying the UKIDSS database 
and removing corrections applied to the aperture photometry in 



SYNMAGs: Catalog-Level Matched Photometry 



9 




0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 

BOSS spec-z BOSS spec-z 



Fig. 4. — For the same sample of BOSS galaxies, the quality of photometric redshifts derived from FWHM-homogenized ugrizY J H K 
GAMA photometry (left panel) compared to the synthetic aperture photometry method using the same filter set (right panel). The SYNMAG 
approach performs better. 



chose to synthesize 2'.'8 radius UKIDSS aperture pho- 
tometry. While we primarily use UKIDSS DR8 catalogs, 
the GAMA photometry was based on DR4, so we note 
that we find no differences when using UKIDSS DR4 cat- 
alogs instead. The SYNMAGs model the varying PSF in 
each of the YJHK bands as a Gaussian with the lo- 
cal FWHM taken from the LAS catalogs (note that a 
more sophisticated PSF treatment is possible). The typ- 
ical UKIDSS FWHM is O'.'S, which is smaller than the 
I'.'l seeing SDSS imaging that defines the intrinsic pro- 
files needed for SYNMAG photorn etry. Extinction correc- 
tions from iSchlegel et al.l ()1998D are applied using val- 
ues recorded in both the SDSS and UKIDSS catalogs. 
Note that the Kron-aperture-matched photometry from 
GAMA DRl is also extinction corrected. 

The restframe (r — J) color-magnitude diagrams for 
both the GAMA photometry (left panel) and SYNMAG 
photometry (right panel) are shown in Figure [31 The 
choice of (r— J) provides a test of the photometric match- 
ing between the SDSS and UKIDSS data sets. To provide 
a reference, the contours indicate the location of much 
less luminous COSMOS galaxieJH with similar redshifts 
as the CMASS sample. We note an offset that may arise 
from the combination of zeropoin t calibration uncer tain- 
ties in the COSMOS data set ^see lCapak et al.ll2007D and 
the fact that the very massive BOSS galaxies are not well 
represented in the small volume of the COSMOS field. 
Because photometric scatter is likely to dominate over 
intrinsic scatter, the observed scatter in these diagrams 
gives some indication of the performance of the photom- 
etry matching. The SYNMAG photometry appears to do 
very well, providing a smaller scatter (cr — 0.32) com- 
pared to the GAMA photometry {a = 0.46). We note 
however that GAMA DRl photometry may have been 
affected by a bug in the PSF convolution (S. Driver, pri- 
vate communication). Our own tests suggest its impact 
on the GAMA colors may not be significant, but future 
releases of GAMA photometry may show improvement. 

the WFCAM Science Archive (WSA). These will be described in 
detail in Bundy et al. (in preparation). 

1^ The COSMOS data set is de scribed in IBundv et~all pOlOl '). 
with restframe colors calculated bv lllbert et al.l l|2010l ). 



3.3. Photometric Redshifts 

The moderate scatter in the SYNMAG SDSS-UKIDSS 
colors of galaxies on the red sequence suggests that syn- 
thetic aperture photometry performs well, but the scatter 
alone does not tell us whether the colors are "correct." 
We therefore perform a second test by comparing the 
performance of photometric redshifts derived using both 
the GAMA and SYNMAG approach. Here the goal of 
accurate photo-zs is more clearly defined, although opti- 
mizing photometry based on the performance of photo- 
zs carries implicit assumptions that should be carefully 
considered (Section |4|). 

We again turn to BOSS galaxies with spectroscopic 
redshifts as a benchmark for evaluating the photo- z per- 
formance. We utilize all BOSS targets in this case (in- 
cluding the LOWZ sample) which allows us to probe red- 
shifts from 0.1 to 0.8, thereby testing the ability of multi- 
ple SDSS-UKIDSS filter combinations to correctly mea- 
sure restframe spectral features. We apply the vl.99.3 
Bayes ian Photometr ic Redshift (BPZ) code (Benite3 
2000; iCoe et al.ll2006l ) to both the GAMA matched pho- 
tometry and the SYNMAG photometry. BPZ includes a 
procedure for estimating photometric zeropoint offsets 
from the mismatch between the best-fit templates and 
the observed photometry of galaxies where the redshift 
is known spectroscopically. The resulting suggested off- 
sets were ~0.02 mag in the u through J bands, but 0.1 
mag in the -ff-band, and 0.18 mag in the JsT-band, re- 
flecting the poorer quality of the BPZ templates in the 
near-IR. We apply these offsets to both sets of photom- 
etry, although their impact is small. 

The comparison of derived photometric redshifts com- 
pared to the 2989 BOSS-measured spectroscopic red- 
shifts is shown in Figure 01 We emphasize that the goal 
is not to obtain the best photo-zs possible for this sam- 
ple, but to compare the resulting photo-z quality from 
two different photometry techniques. When excluding 
outliers defined as |zspcc — -Zphotl > 0.2, we find that 
the FWHM-homogenization approach used by GAMA 
achieves a dz/(l -I- z) of 0.06 while the synmags deliver 
an improved scatter of 0.05. Indeed, the decreased scat- 
ter in the SYNMAG photo-zs is evident by eye in Figure [H 
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The SYNMAGs outlier fraction is slightly larger, however 
(1.7% instead of 1.2%). 

This result was anticipated by Figure |3l the tighter col- 
ors do appear to yield better photo-zs. This is highly en- 
couraging for the utility of the SYNMAG approach. Even 
if future improvements to the GAMA photometry de- 
liver better BPZ performance, the SYNMAGs require far 
less work and computation time. The photometry com- 
putation for the comparison set of ^3000 galaxies here 
takes less than one minute on a modern desktop (8-core 
CPUs with 10 GB memory). Combined with the red- 
sequence test, these results indicate that the scatter in 
SYNMAG colors may outperform PSF homogenization. 
The photo- z test also tells us that SYNMAGs have not 
introduced strong systematic offsets (greater than 0.05 
mag) between filters because these would decrease the 
photo-z precision. The more general question of whether 
SYNMAGs or any matching technique introduces biases 
in the observed colors requires defining what is meant by 
"correct" colors, against which biases could be measured. 
As we discuss in Section |4l no single definition is appro- 
priate and this argues for adopting different approaches 
tied to specific scientific goals. 

3.4. Why SYNMAGs Perform Well 

A comparison of the red-sequence colors (Figure [3]) 
and photometric redshifts (Figure |4]) suggests that sets 
of photometry constructed with SYNMAGs perform bet- 
ter than FWHM-homogenization matched-aperture pho- 
tometry. Putting aside the possibility that the GAMA 
matched-aperture photometry may improve in future re- 
leases, it is worth asking why SYNMAGs may perform bet- 
ter in these tests. Two explanations are apparent. To be- 
gin with, SYNMAGs do not require a degradation of the 
image quality to the lowest-common-denominator PSF. 
For many data sets, this is typically FWHM=2", which 
is on the same order as both the size of distant galaxies 
and the aperture size over which their flux is measured. 
The benefit of not degrading the PSF is demonstrated 
by the improved scatter of the SYNMAG (r — J) colors. 

The second explanation is that we have constructed 
a set of photometry that takes advantage of profile- 
fitting photometry when available. The photo-zs pre- 
sented above are based on a mix of SDSS ModelMag 
colors from profile-fitting in the optical and SYNMAG col- 
ors defined by aperture photometry from UKIDSS in 
the near-IR. We can examine the value of the profile- 
fitting colors alone by ignoring the near-IR SYNMAGs for 
the moment and comparing colors and photo-zs between 
GAMA and SDSS in the optical only. Figure [S] presents 
the restframe {u — r) color-magnitude diagrams derived 
for CMASS galaxies from GAMA (left panel) compared 
to the SDSS profile-fitting ModelMags (right panel), and 
Figure [H] shows the same test with ugriz photo-zs. Both 
show significantly reduced scatter from the use of Mod- 
elMagf[3- The SYNMAGs themselves cannot claim credit 
for the reduced scatter in this case, but it is clearly im- 
portant to make use of profile-fitting photometry if it is 

It is perhaps surprising that the loss of near-IR photometry 
does not degrade the photo-zs. In fact, for the ModelMags, the 
photo-zs significantly improve without the near-IR. This is in large 
part a result of uncertainties in restframe near-IR templates, a 
widespread problem that has motivated template error weighting 
in recent photo-2 codes (e.g.. [Erammer et al.ll2008 ). 



available. It is likely, for example, that a mix of Model- 
Mags in the optical and FWHM-homogenized colors in 
the near-IR would yield tighter colors and photo-zs than 
FWHM-homogenized colors in all bands. 

4. CONCLUSIONS 

We have presented a technique for using catalog-level 
data products to match the photometry of disparate 
imaging data sets. Our SYNMAG approach uses multi- 
Gaussian expansions to represent both intrinsic profiles 
and the PSF, allowing for extremely fast convolutions. 
The result is a tool that can be used to quickly ob- 
tain matched photometry for large data sets without re- 
analyzing pixel data. 

We find that SYNMAGs yield colors and derived photo- 
zs with scatter that is as good if not better than FWHM- 
homogenized aperture photometry. We have also shown 
that the resulting scatter is minimized by constructing 
photometric data sets that employ profile-fitting colors 
when available. These tests, however, assume that the 
"best" color estimators minimize the scatter. This may 
not be the case, especially given that most galaxies are 
2-coniponent bulge-|-disk systems, often with different 
colors for each component. In fact, the choice of aper- 
ture size or assumed profile shape should depend on 
the goal for which the photometry is being analyzed. 
Such "application-oriented photometry" would suggest, 
for example, that for photo-zs, fitting profiles that em- 
phasizes the typically redder light from the bulge would 
strengthen the 4000A break in the derived SED and yield 
tighter photo-z estimates. The same profile shape would 
not be optimal for studying the SEDs of galaxy disks, 
however. 

Applying SYNMAGs to future survey data would bene- 
fit from having some specific measurements recorded in 
publicly available catalogs. Even if profile-fitting pho- 
tometry is performed, this basic information is valuable 
for cross-checks, troubleshooting, and the application of 
SYNMAGs to future overlapping data. First, it would be 
useful to have aperture photometry performed in a series 
of circular apertures that range from roughly one-half of 
the typical PSF FWHM to i? = 6-7" (integer values like 
radii of 2", 3", 4", etc., are common). The R = 6" cor- 
responds to ^11 kpc//i7o at z = 0.1, sufficient for galaxy 
photometry at z > 0.1. Several larger apertures up to 
R — 20" would provide adequate sampling for most of 
the largest galaxies in SDSS (z < 0.1). No corrections 
or adjustments (e.g., PSF corrections) should be applied 
to the aperture fluxes beyond background subtraction. 
Reliable uncertainties on each flux measurement are also 
valuable. Special care is needed for smaller apertures 
that subtend only a few pixels, otherwise it may be better 
to remove small aperture measurements from the cata- 
logs. It is also useful to consider the definition of null val- 
ues for measured magnitudes. The user should be able to 
distinguish between the case when no flux was detected 
(e.g., TO — 99) and when imaging for a source in a par- 
ticular band is missing (e.g., to — —99), either because 
it was not observed or was somehow compromised. 

SYNMAGs also require information about the PSF. Ide- 
ally, each catalog entry would provide a way to deter- 
mine the PSF shape at the location of the measured ob- 
ject given the conditions under which it was observed. 
Often a 2-component Gaussian description is sufficient 



SYNMAGs: Catalog-Level Matched Photometry 



11 



I 2: 



: 

GAMA 
\ a = 0.79 








• . % 

• .••'-..•V.''« 


^iVvfi r. 










»: .V-.'. '• 




' ■ 





-21 



-22 



-23 



-24 



4 


1 ■ 

ModelMags 
O — u.o^ 


1 1 1 1 1. 


3 




• * * " 

. : ^ T .' • 

• • • . 


1 2 


• 














1 




"... ■" 





I ■ 





-21 



-22 



-23 



-24 



Fig. 5. — Restframe (n — r) color-magnitude diagrams of (mostly) red-sequence CMASS galaxies from the BOSS survey. As with Figure[3l 
the observed scatter for the GAMA FWHM-homogenized matched-aperture photometry is larger than for SDSS ModelMags, demonstrating 
the utility of profile-fitting in this context. Contours indicate the location of much less luminous COSMOS galaxies in the same redshift 
range. 




ahhough the modeled shape could be more sophisticated 
and could be modeled as a function of both detector posi- 
tion and observing conditions or exposure number. The 
position angle of off-axis or non-axisymmetric compo- 
nents also needs to be recorded. 
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APPENDIX 

COOKBOOK: USING THE SYNMAG SOFTWARE 

This appendix presents instructions for using the SYNMAG softwar43. We describe the key steps in the process and 
point to the routine synmag_sdss for the specific case of matching to SDSS data. 

Before using the SYNMAG routines, some photometric information from the data sets involved is required. This 
should include parameters from profile fitting in at least one band, aperture photometry results for the "target" data 
set where profile fits are not available, and seeing estimates for this data set. With these information in hand, the 
procedure is as follows: 

1. Decompose the intrinsic profile shape into an MGE. This only needs to be done once for each profile shape that 
is used. For standard de Vaucouleurs and exponential profiles, or their SDSS variants, the Gaussian coefficients 
can be read off of Table [TJ Assuming more generally that the profile is a Sersic, the mgsersic routine can be 
used to determin e the MGE. It utilizes a ID fitting routine in the mge_fit_SECTORS IDL package presented in 
iCappellaril (|2002[ ). As disc ussed in Section 12.4. 1[ however, a least-squares MGE fit to the 2D profile provides a 
better approximation fsee lHogg &: Landl2012D . The PSF of the target catalog can also be decomposed in this 
way if desired. 

2. Compute synthetic photometry. Use synmag to compute the synthetic aperture photometry in the profile-defining 
band for the PSF and aperture size measured in the target catalog. Multiple apertures can be requested and the 
PSF parameters can be specified separately for each source to account for variations across the image or from 
differing observing epochs. 

3. Assemble a matched set of photometry. Repeat the synmag procedure for each filter band in the target catalog 
to account for PSF variations and missing data from one band to the next. Then subtract the resulting colors 
from the total flux measured in the profile-defining band in order to construct a set of matched flux total flux 
measurements (i.e., Wi^tot = Wp_tot — Cpi, see Section [2TT|) . 

The PSF-convolved MGEs can also be returned by synmag. The routine mgestr computes the effective radius and 
axis ratio of any MGE (see Section 12.4. 3p . 
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